Fluxon readout of a superconducting qubit 
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Quantum computing using superconducting 
circuits underwent rapid development in the last 
decade^^. This field has propelled from quan- 
tum manipulation of single two-level systems 
to complex designs employing multiple coupled 
qubits allowing one to execute simple quantum 
algorithms^i^. On the way to a practical quan- 
tum computer, a need for scalable interfaces be- 
tween classical circuits and the quantum coun- 
terparts has arisen. Low-temperature supercon- 
ducting single-flux quantum (SFQ) logic employs 
magnetic fluxons in Josephson transmission lines 
(JTLs) as basic bits for classical computation^'*. 
Here, we report on an experiment implementing 
a direct link between SFQ electronics and a su- 
perconducting qubit. We demonstrate a readout 
of the state of a flux qubit through a frequency 
shift of a single fluxon oscillating in a JTL. The 
energy spectrum of the flux qubit is measured us- 
ing this technique. The demonstrated approach 
may open ways to future full-scale integration of 
solid-state quantum computers with digital SFQ 
electronics. 

A quantum computer requires a set of coupled quan- 
tum bits (qubits) with control gates to manipulate them, 
as well as a classical computer as an interface to the 
quantum counterpart^. As for today, the most promis- 
ing candidates for the role of quantum bits come from 
the field of solid state superconducting qubits - as phase, 
flux or transmon qubits^"— 1^°'^"'^ . Experiments of the 
past decade have shown how to entangle and operate 
these qubits^ii^, and simple quantum algorithms have 
been demonstrated^. Nowadays, a classical interface 
for qubits is an emerging milestone in the development 
of circuits with multiple solid state qubits. An efficient 
control and readout of several quantum bits requires a 
powerful classical computer in order to process the vast 
amount of real-time measurement data from the quan- 
tum counterpart. Researchers often use specific pro- 
grammable electronic boards, e.g., field-programmable 
gate array (FPGA) boards, in order to meet the require- 
ments for high processing speed, flexibility and reason- 
able price. However, in the near future this approach 
will become inefficient due to the design complexity and 
engineering problems of communication between many 
qubits and room temperature readout electronics. 



Superconducting SFQ electronics offers a possibil- 
ity of implementing an extremely fast, up to 500 
GHz clock speed, low-temperature classical computer— 1^. 
SFQ logic can be optimized to be used at millikelvin 
temperatures^"*'^- and be also modified to allow for re- 
versible operatioE^iii. The lack of compact memory el- 
ements in SFQ electronics can be resolved by using mag- 
netic Josephson junctions*®. With the advent of super- 
conducting quantum computing and the requirement to 
process a lot of data at low temperatures, SFQ electron- 
ics seems to be the natural candidate for the role of an 
interface between room temperature electronics and its 
quantum counterpart. In this paper, we experimentally 
demonstrate the key link between a superconducting flux 
qubit and SFQ-based electronics. 
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FIG. 1. An annular Josephson junction coupled to a flux 
qubit. The junction is homogeneously biased by an external 
current source It . The flux qubit is controlled by a dc control 
line current Icl and a microwave signal lac at the frequency 
oive- The purple lines schematically show the flow of the elec- 
trical currents. The fluxon equilibrium oscillation frequency 
i/Q is shifted by an amount Sv due to fluxon scattering on the 
current dipole produced by a flux qubit. Microwave radia- 
tion from the fluxon is picked up by a capacitively coupled 
microstrip antenna and later fed to a cryogenic amplifler (not 
shown) . 



A single fluxon (alias Josephson vortex) in an under- 
damped JTL has properties of a relativistic particle car- 



rying a magnetic flux quantum $0 = h/2e^. The size 
of a vortex can vary from a few to several hundreds of 
microns, depending on the critical current density jc and 
vortex velocity in JTL. By applying a bias current, the 
vortex can be accelerated up to the Swihart velocity cs^ 
which is the characteristic velocity of electromagnetic 
waves in JTL. The dynamical properties of the fluxon 
resemble a classical particle with a well-defined mass and 
velocity. In SFQ technology, fluxons are used to realize 
logic operations by transmitting and detecting them in 
JTLs. 

The idea of using fluxons for detecting the state of 
flux qubits was initially proposed and theoretically ana- 
lyzed in Refi^. A ballistic fluxon moving in a JTL that 
is weakly coupled to a superconducting flux qubit can 
be used to readout the quantum state of the qubit^i^i. 
There is a potential advantage of the fluxon readout 
over the usual dispersive readout of qubits with the use 
of microwave resonators, usually referred to as QED 
readout^. The distribution of the electro-magnetic field 
in the usual bi-coplanar resonator is such that it per- 
turbs the qubit as long as there are photons in the res- 
onator. For the fluxon readout, the back-action on the 
qubit occurs only for a very short time, when the fluxon 
passes by the qubit. This time depends on JTL parame- 
ters with the typical time scales of about a few picosec- 
onds. This short-time interaction improves isolation of 
the qubit from the environment and hence makes it pos- 
sible to increase the qubit coherence. 

For our experiments, we employ a JTL formed by a 
continuous annular Josephson junction (AJJ). JTLs of 
such closed topology (see Fig.[T]) are used as on-chip SFQ 
clocks^'^ and their most significant advantage is the con- 
servation of the total magnetic fiux initially trapped in 
JTL. One can create a fluxon on demand in a flux-free 
JTL by applying a current through a pair of injectors^^. 
To couple a flux qubit to the fluxon inside the AJJ, it is 
necessary to engineer the interaction between two orthog- 
onal magnetic dipoles. To facilitate this interaction, we 
have added a superconducting coupling loop embracing 
the flux qubit, as shown in Fig. [1] The current induced 
in the coupling loop attached to the AJJ is proportional 
to the persistent current in the flirx qubit. Thus, the per- 
sistent current Ip in the qubit manifests itself in the AJJ 
as a current dipole with an amplitude /^ = kIp/(jcXjW) 
on top of the homogeneous background of bias current 
7 — h/Ic- Here, the coefficient k reflects the inductive 
coupling between the qubit and the coupling loop as well 
as the fluxon differential resistance, W is the width of 
the JTL, Xj is the Josephson penetration depth, jc is 
the critical current density, /{, is the bias current, and 
Ic denotes the critical current of the AJJ. A theoreti- 
cal description of interaction between a single Josephson 
vortex and a current dipole in the AJJ can be given by 
the perturbed sine-Gordon equation (PSGE)?^ 



with the corresponding boundary conditions 
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where n is the number of trapped fluxons, a = ujp/uJc 
is the damping parameter, I = L/Xj is the normal- 
ized junction circumference, d = D /Xj is the normalized 
dipole length, ujp = ^j2eI^JhC is the plasma frequency, 
ujc — 2eIcRN /fi is the characteristic frequency, C is the 
AJJ capacitance, Rm is its normal state resistance. 

Under the influence of the bias current /;,, the fluxon 
rotates in the AJJ. This rotation produces both dc and 
ac voltage components, with the magnitude of the for- 
mer proportional to the frequency of the latter. In our 
detection scheme, we directly measure the ac voltage in- 
duced by the fiuxon as microwave radiation using a low- 
noise cryogenic microwave amplifier^^. This approach to 
measure the fluxon radiation frequency provides a much 
greater precision than dc voltage measurements. 

The studied JTL was fabricated by using photolithog- 
raphy and the standard Nb/A10a;/Nb trilayer process 
with critical current density jc — 800 A/cm^. The es- 
timated Josephson penetration depth is Xj ~ 13 /^m, the 
Josephson plasma frequency is ujp/2iT ~ 154 GHz, and 
the estimated damping parameter at the operating tem- 
perature is estimated as a ~ 0.02. The circumference 
of the junction L = 880 /im determines the maximum 
frequency of the radiation at about 15 GHz, correspond- 
ing to the fluxon moving with the Swihart velocity cs- 
The width of the AJJ was W — 5/im, the dipole length 
D ~ 35 /im and its fluxon-free critical current Ic — 35 
mA. The flux qubit was deposited using the standard 
aluminum shadow evaporation proces^iS on top of the 
pre- fabricated niobium structures. Designed fabrication 
parameters for the Josephson junctions in the fiux qubit 
loop were the following: critical current Ic = 500 nA, 
alpha factor aq — 0.63, ratio of charging and Josephson 
energies Ec/Ej = 0.0034. 

To experimentally test the qubit readout scheme, the 
temperature was set at T = 70 mK, well below the super- 
conducting transition temperature Tc of aluminum form- 
ing the qubit. The AJJ was biased at a fixed current 
/{,. Then, we varied the current through the control line 
IcL in order to change the magnetic flux through the 
flux qubit. Due to periodic variations of the persistent 
current Ip in the qubit loop versus the control line cur- 
rent IcL, the current dipole strength /j, is also changing 
periodically. The persistent current of the flux qubit in 
the ground state was calculated numerically^ and is de- 
picted by the black solid line in Fig. [21a. 

When a ballistic fiuxon scatters on a positive current 
dipole, it first gets accelerated and then decelerated by 
the dipole poles. In the case of absence of damping and 
bias current, the change of its circulation frequency is 
determined only by the dipole polarity and the dipole 
amplitude. In the presence of finite damping a and ho- 
mogeneous bias current 7, the total propagation time be- 
comes dependent on the complex interplay between bias 
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FIG. 2. a) The theoretical persistent current Ip of the 

ground state of the flux qubit versus magnetic frustration 
(black line) calculated for our flux qubit parameters. The 
red line shows the expected fluxon frequency shift in kHz for 
7 — 0.4 and n = 0.05. b) The experimentally measured mod- 
ulation of the fluxon oscillation frequency due to the coupling 
to the flux qubit. Black dots show the measured mean fre- 
quency of fluxon oscillations. Every point consists of 100 av- 
erages; bias is 7 = 0.39. The red line shows the corresponding 
fit. 



current, current dipole strength, and damping. Our pre- 
vious numerical simulation of Eq. ([T]) and perturbation 
analysis^^ have shown that, for fluxon velocities close to 
cs, the shift of fluxon oscillation frequency almost does 
not depend on the sign of the dipole fi and depends only 
on the absolute amplitude of the dipole /i and the bias 
current 7, looking like it is shown by the red line in 
Fig. [2ja (where the black curve Ip was taken as the cur- 
rent dipole strength ~ /u). As the dipole length d ~ 3 
is much larger than the characteristic size of the fluxon 
at relativistic velocities, the contributions of the sepa- 
rate dipole poles to Sv is additive and is not dependent 
on their order. For smaller bias currents (small fluxon 
velocities) and smaller dipole lengths, when d is compa- 
rable with the fluxon size, Su becomes dependent on the 
dipole sign. 

The back-action of the propagating fluxon on the flux 
qubit can be estimated from the magnetic flux <i>f,a which 
the fluxon excites in the qubit loop, depending on the 
geometrical inductances of the loop and the qubit. For 
our test sample, the estimated maximum back-action flux 
^ba is about 35 m$o- Considering that the fluxon inter- 
acts with the qubit on the time scale of r ~ 3 ps, we 



can calculate the qubit phase shift S(f) = L 



r Ei2(t)~Et 



^dt 



FIG. 3. a) The theoretically expected persistent current Ip 
for the ground state of the flux qubit versus magnetic flux 
(black line). The spikes indicate the presumed transitions 
to the first excited state at flux bias points corresponding to 
E\2 = 14 GHz flux qubit energy splitting. The green line 
shows the expected persistent current Ip for the first excited 
state, b) Modulation of the fiuxon's oscillation frequency 
due to the coupling to the qubit, measured in the power do- 
main. Measured power P relates to the power at the fixed 
frequency offset -1-50 kHz from the fluxon mean oscillation 
frequency Vf. Every point consists of 10 averages with video 
filter bandwidth of 10 Hz. The red line depicts the corre- 
sponding theory fit. 



induced by the fluxon. We find that 6(t> ^ 2- 10^"^ • 2-k per 
fluxon revolution in the AJJ. It corresponds to the de- 
phasing time due to the readout back-action of T^a '^ 35 
ns which is consistent with the observed linewidth of the 
flux qubit spectra. This short coherence time may be 
considered as a strong disadvantage of the fluxon read- 
out. However, simply by decreasing the critical current 
density to jc ~ 10 A/cm^ and, at the same time, by in- 
creasing the length of the coupling loop, one should be 
able to decrease the back-action and increase T^a by, at 
least, 2-3 orders of magnitude keeping the current dipole 
amplitude /i at the same level. 

Figure [21b shows the experimental shift of the fluxon 
oscillation frequency due to coupling to the flux qubit 
at the fixed bias current 7 = 0.39. We measured a 
Lorentzian radiation peak of fluxon oscillations in the 
AJJ for fixed values of Icl and then determined the mean 
frequency vj for which the radiation power was maximal. 
To fit our experimental data to the theory, we take into 
account the parasitic coupling between the control line 
current and the fluxon leading to an additional linear 
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FIG. 4. Fluxon response on the microwave excitation applied 
to the flux qubit. Color scale corresponds to the detected 
power. Measured power P relates to the power at the fixed 
frequency offset +50 kHz from the fluxon mean oscillation 
frequency Vf. Every point consists of 10 averages with video 
filter bandwidth of 1 Hz. The flux qubit spectrum can be 
clearly seen as the blue curved trace. The black dashed line 
depicts the corresponding fit of the flux qubit. 



shift of Vf versus Icl ■ The corresponding fit is presented 
by the red line in Fig. [2lb showing fairly good correla- 
tion between theory and experiment. Smaller irregular 
peaks visible in Fig. [21b are low-frequency fluctuations, 
most probably arising from trapped Abrikosov vortices in 
the superconducting leads of the A J J. Presumably, these 
Abrikosov vortices are also responsible for the parasitic 
flux offset in Fig. [2]b as this offset varies for different 
cooldowns. 

To speed up the measurement process and decrease 
noise, we switched to zero-frequency span and measured 
the power at the fixed frequency vj + 50 kHz. Fig- 
ure [31b shows the results of such measurements around 
$ ~ $o/2- A sharp dip is clearly visible at the $o/2 flux 
as well as two other smaller symmetric satellites around 
it. These satellites occur when the energy splitting Eqi 
between the ground and excited states of the qubit coin- 
cides with the fluxon oscillation frequency Vf . This leads 
to the change of persistent current in the qubit for two 
values of flux bias $ shown by the black line in Fig. [31a. 
The theoretical fit of the data for the fluxon resonance 
frequency Vf = 14.031 GHz is presented by the red line 
in Fig. [3lb. At resonance, pumping at the fluxon fre- 
quency should lead to Rabi oscillations in the flux qubit 



between its ground and excited states and, therefore, the 
measured signal reflects a mixture of the ground and the 
first excited states. This operation is similar to resonant 
interaction between qubit and resonator in a cavity QED 
setup. 

As the last step, we performed microwave spectroscopy 
of the qubit. We applied an additional continuous exci- 
tation tone from an external microwave source at the 
frequency Ve and swept the control line current between 
two resonant dips (see Fig.[3lb). The signal of the fluxon 
readout for every flux bias point without microwaves was 
subtracted from the actual response with microwaves to 
get rid of an unwanted background slope. The result- 
ing color plot of P{h'e,IcL) is shown in Fig. [H One 
can clearly recognize the hyperbola of the flux qubit 
spectrum^^ as a white-blue curved line between 2 and 
10 GHz, with the minimal energy splitting A ~ 2.7 
GHz. We can very well flt the measured spectrum by 
the theoretical curve (shown by the black dashed line 
in Fig. [3]) for the following parameters: critical current 
Ic — 320 nA, alpha factor aq ~ 0.58, ratio of charging 
and Josephson energies Ec/Ej = 0.0034. It is worth 
noting that the fluxon readout in this experiment was 
operated in the classical rather than quantum regime. 
Employing fluxons to interact with or act as qubits in 
the quantum regime, as discussed theoretically^S- and also 
detected experimentally22i2£^ would be very exciting and 
challenging in the future. 

In conclusion, we have demonstrated the coupling of 
an oscillating single fluxon in an annular Josephson junc- 
tion to a flux qubit. The persistent current in the flux 
qubit has been detected as a shift of the fluxon oscilla- 
tion frequency. Resolution of the measurement scheme 
is high enough to measure transitions between quantum 
states as shown by the acquired energy spectrum of the 
flux qubit. The sensitivity of this readout can be signifi- 
cantly improved in various ways, for instance, by lowering 
the critical current density jc or decreasing the width W 
of the long Josephson junction, as these both increase 
the current dipole strength and therefore the fluxon fre- 
quency shift. Our results prove the possibility of detect- 
ing quantum states by fluxon readout and, thus, open 
the way to applying SFQ technique in the fleld of super- 
conducting quantum computing. 
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